ASTROPHYSICAL JOURNAL, IN PRESS (ACCEPTED 2013 MAR 04). 
Preprint typeset using KTgX style emulateapj v. 08/13/06 



SN 2009JS AT THE CROSSROADS BETWEEN NORMAL AND SUBLUMINOUS TYPE IIP SUPERNOVAE: OPTICAL 

AND MID-INFRARED EVOLUTION 

P. Gandhi 1 - 2 , M. Yamanaka 3 A5 , M. Tanaka 6 , T. Nozawa 7 , K.S. Kawabata 5 , I. Saviane 8 , K. Maeda 7 , T.J. Moriya 7 - 910 , T. 

Hattori", M. Sasada 12 , R. Itoh 3 

Astrophysical Journal, in press (Accepted 2013 Mar 04). 

ABSTRACT 

We present a study of SN 2009js in NGC 918. Multi-band Kanata optical photometry covering the first ~ 120 
days show the source to be a Type IIP SN. Reddening is dominated by that due to our Galaxy. One-year-post- 
explosion photometry with the NTT, and a Subaru optical spectrum 16 days post-discovery, both imply a good 
match with the well-studied subluminous SN 2005cs. The plateau phase luminosity of SN 2009js and its plateau 
duration are more similar to the intermediate luminosity IIP SN 2008in. Thus, SN 2009js shares characteristics 
with both subluminous and intermediate luminosity SNe. Its radioactive tail luminosity lies between SN 2005cs 
and SN 2008in, whereas its quasi-bolometric luminosity decline from peak to plateau (quantified by a newly- 
defined parameter Alog£ measuring adiabatic cooling following shock breakout) is much smaller than both the 
others. We estimate the ejected mass of 56 Ni to be low (^0.007 M Q ). The SN explosion energy appears to have 
been small, similar to SN 2005cs. SN 2009js is the first subluminous SN IIP to be studied in the mid-infrared. 
It was serendipitously caught by Spitzer at very early times. In addition, it was detected by WISE 105 days 
later with a significant 4.6 /im flux excess above the photosphere. The infrared excess luminosity relative to the 
photosphere is clearly smaller than that of SN 2004dj extensively studied in the mid-infrared. The excess may 
be tentatively assigned to heated dust with mass ^3 x 10 -5 M , or to CO fundamental emission as a precursor 
to dust formation. 

Subject headings: supernovae — individual: (SN 2009js, SN 2008in, SN 2005cs) 



1. INTRODUCTION 

Core-collapse supernovae (SNe) represent the deaths of 
massive stars through runaway implosions of cores, followed 
by explosive expulsion of their outer layers. They exhibit a 
diverse range of properties in terms of their spectral and pho- 
tometric evolution (e.g. Filippenko 1997[ ). Observational 
follow-up within the past few years has also uncovered an 
incredibly broad range of emitted luminosities: from super- 
luminous SNe ( |Gal-Yam]2 012) with luminos ities greater than 



7x10 erg s 1 , to 'subluminous' events (Pastorello et al. 



2004 1 that are factors of 10-100 times less powerful at peak. 
Understanding the origin of this diversity is an important goal 
of present research. 
The Type II class of SNe events display spectroscopic sig- 
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natures of Hydrogen, and appear in at least three sub-classes. 
The IIP sources are named after a stage of flat ('Plateau') light 
curves showing constant flux for several months. IIL sources, 
on the other hand, show a monotonic and linear decline of 
their post-explosion brightness without any plateau. Lastly, 
the Iln sub-class shows spectral features that are much nar- 
rower than in the other sub-classes, thought to be a result 
of strong interaction between SN ejecta and the circumstel- 
lar medium. Other sub-classes include the lib (with weak 
hydrogen spectral features that soon disappear), and ambigu- 
ous events collected under the generic H-pec, or peculiar, sub- 
class. 

The IIP population is important because it represents the 
most numerous sub-class, constituting about 75% of all Type 
II SNe ( |Li et al.|201 T]|Smith et al.|201 1) . This sub-class itself 
encompasses a wide range of properties (Hamuy |2003| l, and 
there is much debate regarding the progenitors of these events 
with prese nt constraints f avoring masses of between 8~16 
M Q (e.g. Smartt et al.|2009 for a review). The subluminous 
SNe may occupy the lower mass end of this range ( |Fraser| 
et al. |201 1 1, though comparatively few firm constraints ex- 



ist. Their low luminosities have been attributed either to low 
mass progenitors, or to 'fall back' events in which a massive 
progenitor gives rise to a black hole core and swallows a large 
fraction of the stellar material which would otherwise have es- 



caped, thus limiting the observed luminosity (e.g. Zampieri 
et al.|1998| Moriya et al.|20 10). Only about a dozen sublumi- 



nous IIPs have been studied in any detail (e.g. |Turatto et al. 
[T998l |Pastorello et al.|2004"llSpiro & Pastorello|2009[ IFraser 
|et al.|[201 l| l, and understanding their nature will require ob- 
servations of more sources to enlarge the present sample and 
also to bridge the 'ga p' between the subluminous and normal 
IIP population (e.g. Roy et al.|201 1) >. 

In this paper, we present optical and mid-infrared obser- 
vations of SN 2009js in the galaxy NGC 918. As we show in 
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this work, the source was a Type IIP event, which shares many 
characteristics with two other SNe: the well-studied sublumi- 
nous SN 2005cs, and the intermediate luminosity SN 2008in. 
We report optical photometric evolution over the first year of 
SN 2009js, in addition to optical spectroscopy early in the 
plateau phase. From the multi-band light curves, we can con- 
strain the explosion epoch, and measure the plateau length 
and source luminosity, thus allowing a computation of the 
line-of-sight reddening and explosion energy. Mid-infrared 
(MIR) photometric data with the Spitzer and WISE missions 
are also reported, and allow tentative constraints on the dust 
surrounding SN 2009js. As far as we are aware, this is the 
first mid-infrared study of a SN which shares characteristics 
of subluminous events. This work is thus a valuable addition 
to the low luminosity class of events. 



The discovery of SN 2009js was first reported by Nakano 
2009), and also independently by Silverman et al. 



et al. 

p009l > dating from 2009 Oct 11.44 and Oct 11.689. The 
SN was located at J2000 coordinates of RA=02h25m48.3s 
Dec=+18°29'26", about 35"5 W and 2077 S of the nucleus 
of NGC 918. Optical spectra obtained the following day with 
the 3-m Shane reflector (+ Kast) at Lick Observatory indi- 
cated the object to be very similar to the Type IIP SN 2005cs 
about two days post-maximum brightness. An Ha absorption 
blueshift of 7200 km s _1 intrinsic to the SN was measured in 
this spectrum. 
We assume 

=0.005 



a distance to the host galaxy NGC 918 of 
21.7±1.8 Mpc O 



e.g. 



Martinez-Garcia et al.l 2009 



and references therein). This galaxy is of morphological class 
Sc 13 oriented at a position angle of 158°. It has hosted two 
SNe in the past three years: SN 201 lek (a Type la reported 
by |Nakano et al. 201 l| l, and SN 2009js which is the subject of 
this paper. 

2. OBSERVATIONS 
2.1. Optical imaging and spectroscopy 
Three sets of optical imaging data are used in our work: 

1. Early -time photometry, including pre-discovery limits 
and d iscovery data (Nakano et al. 2009; Silverman et al. 
2009). These were obtained in filterless imaging mode. 
Throughout this paper, we assume a discovery epoch of 
2009 Oct 11.44 ( |Nakano et aLl[2009l [Silverman et al.| 
|2009l >. 

2. Plateau-phase BVRqIq follow-up of the field of 
SN 2009js earned out with the HOWPol inst rument 
(Hiroshima One-shot Wide-field Polarimeter; |Kawa-| 

. |2008] > at the Kanata telescope, starting on 



bate et al. 



2009 Oct 14.8 at an epoch of +3.4 days (epochs are 
quoted post-discovery throughout) and continuing for 
about 4 months. 

3. A single late-time epoch in the VRI filters is pub- 
licly available from observations made at the ESO Faint 
Object Spectrograph and Camera (EFOSC; Buzzoni 
| e t al. 1984) mounted on the New Technology Telescope 
(NTT). These data date from 2010 Oct 6, just under one 
year post-explosion. 

All observation details (including logs and analysis) may 
be found in the Appendix. The photometric measurements 
are listed in Table I All 



In addition to photometry, an optical spectrum of SN 2009js 
was obtained by some of us on UT 2009-10-27 (MJD = 
55131.6), about 16 days post-discovery, with the Faint Object 
Camera and Spectrograph (FOCAS; Kashik awa et al.||2002| l 
mounted on the Subaru telescope. This provides important 
information on the SN type and ejecta kinematics. Observa- 
tional details may again be found in the Appendix. 

2.2. Mid-infrared imaging 

Although the source had been monitored in the optical since 
2009, we identified it as being of particular interest in 201 1 
when we noticed its detection in the mid-infrared serendipi- 
tously by two missions. MIR detections of SNe (especially of 
the low luminosity class that SN 2009js belongs to) are rela- 
tively scarce. 

Spitzer observed the host galaxy NGC 918 for unrelated 
science 14 on two occasions: about one month pre-explosion 
and also two days post-discovery. SN 2009js shows a bright 
detection in the latter set of images in both observed bands 
centered on 3.6 and 4.5 pm (hereafter, channels IRAC1 and 
IRAC2, respectively) of the InfraRed Array Camera (IRAC; 
|Fazio et al.|2004l i. 

On longer timescales, the Wide-field Infrar ed Survey Ex- 
plorer (WISE) satellite (Wright et~aL] |2010| l detected the 
source in its Wl and W2 bands centered on wavelengths of 
«3.4 and 4.6 pm respectively, at an epoch of about +107 days. 
These data are publicly-available as part of the mission all-sky 
survey release. Later non-detections at epochs of about +295 
and +470 days are also available. 

The MIR observation log may be found in the Appendix 
(Table |Cl| l along with a full description of the data analy- 
sis and photometry procedures. MIR flux measurements are 
listed in Table |C2| Spitzer and WISE images are shown in 
Figs. [T] and [2] respectively. 

3. RESULTS 
3.1. Reddening due to the Galaxy 

The optical data were dereddened assuming a standard in- 
terstellar extinction law with Ry=3.l ( Fitzpatrick| 1999] ). The 
Galactic extinction along this line-of-sight is determined to be 
A v =0.95 from the work of |Schlafly & FinkbeinerlpJTT) . A 
10% uncertainty is assumed on this value. 

In the mid-infrared, several more recent determinations 
of the dereddening law have been published, and we used 
an average extinction law computed from these, resulting 
in A 3A um /A v =0.0 1 and A AS Atm M v =0.05 (See § 2.2.3 of 
Gan dhi et al.|201 l] for details and references). 



3.2. Reddening local to the host galaxy 

Extinction local to the source was estimated by using the 
relationship between host galaxy reddening and plateau V— 
I color identified by Olivares et al. (2010). The V—I color is 
that measured at -30 days from the mid-point of the transition 
phase (tpx). The optical light curves and tpx measurement is 
described in the following section, ipx-30 corresponds to day 
80, at which time the V-I color of SN 2009js (corr ected for 
Galactic reddening) is 0.73±0. 13. Using Eq. 7 of Olivares 
leTaLl ( 12010) returns A v (host)=0.18±0.38. 

3.3. Optical light curves 



13 http://leda.univ-lyonl.fr 



14 Cycle 9 PROGID 61060, P.I. K. Sheth. 
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FIG. 1. — Spitzer IRAC1 and IRAC2 pre-explosion images from 2009 Sep 06 at an epoch of 34.8 days before discovery (left column), and 2.2 days post- 
discovery on 2009 Oct 13 (right column). The position of SN 2009js is marked by the green circle. The size of each image is about 2.4'x2', with North being 
up and East to the left. The SN is clearly visible just two days after explosion in the column on the right in both bands. 



The change in brightness between the early-time pre- 
discovery and the discovery epoch photometry implies a 
brightening by at least 1.7 mags in less than 11 days (see 
Table |A1 [ >. Fig. [3] shows the long-term source optical light 
curves. The early-time photometric detections from published 
filterless photometry match our first Kanata V band measure- 
ments closely. The long steady light curves in several (redder) 
Kanata filters are characteristic of Type IIP supernovae. The 
late-time NTT magnitudes probe the nebular phase and show 
a dimming with respect to the peak brightness measured by 
Kanata of about 6.2, 5.2 and 4.8 mags in V, R and /, respec- 
tively. 

In order to characterize the light curve evolutionary phases, 
we fitted the latter part of the / band light curve (beyond an 
epoch of +60 days) using a two component Fermi-Dirac func- 
tion plus lin ear function following the procedure of |Qlivares| 
et al.lpOTO} see their Eqs. 1 and 2). The / band was chosen 
because the comparison with the light curves of SN 2005cs 
in Fig. [3] shows that this band probably covers the bulk of the 
transition phase betw een t he plateau and tail (this will be fur- 
thered discussed in § 4.4 1. Although our Kanata data do not 



cover the beginning of the tail, the NTT photometry allows 
an approximate (i.e. two-point) estimation of the evolution 
during this phase, under the assumption that the Kanata data 
cover the bulk of the transition. The best-fit parameters for 
the transition phase are found to be 



a = 2.24± 1.04 mag 
tpT = 111.3 ± 3.6 days 
w = 4.5 ± 2.4 days 
Po = 0.009 ± 0.004 mag day" 



(1) 
(2) 
(3) 
(4) 



m = 18.46 ± 0.93 mag (5) 

X 2 /dof = 1.7/5 (6) 

where a represents the plateau brightness above the linear 
phase (in mag), ipx denotes the mid-point of the transition 
phase in days (from discovery) and wo is the width of this 
phase (in days). The slope and intercept of the linear radioac- 
tive nebular decline are po and mo, respectively. The good- 
ness of fit is quantified by x 2 /dof, with dof being the degrees 
of freedom. 

3.4. Mid-Infrared photometry 

The MIR detection of SN 2009js was first noticed by us dur- 
ing an examination-by-eye of the WISE all-sky preliminary 
data release fields (first-pass processing prior to the all-sky 
release) of supernovae to have occurred during, or before the 
beginning of, the WISE survey. Thereafter, we checked and 
found the source to be detected with Spitzer as well. 

The Spitzer detection is very significant at both the observed 
wavelengths of 3.6 and 4.5 /im two days post-discovery. As 
compared to the pre-explosion images obtained about 37 days 
prior to this, the flux at the position of the SN increased by 
factors of a t least 25 and 1 1 in the two channels, respectively 
(Table C2 1. This is undoubtedly a strong lower limit to the 



change in flux of the SN itself. 

The WISE data show the source to be clearly present in the 
Atlas images taken about 100 days post-discovery. Table |C2| 
implies a decline in flux by factors of 1.8 and 1.6 (or magni- 
tudes of 0.6 and 0.5) in Wl and W2 over this period, as com- 
pared to Spitzer IRAC1 and IRAC2, respectively. The WISE 
detection is highlighted in Fig.[2]by subtracting smoothed and 
resampled IRAC1 and IRAC2 pre-explosion images from the 
Atlas Wl and W2 images, respectively. Proximity to the host 
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FIG. 2.— (Top and Middle) WISE 3.4 and 3.6 /xm (Wl and W2) Atlas L3 images of NGC 918 taken around 2010 Jan 26 (epoch ~107 days post-discovery) are 
shown in the left column. In the right column, Spitzer pre-explosion IRAC1 and IRAC2 images have been sub tract ed, after scaling and PSF matching. The SN is 
clearly visible. (Bottom) Coadded images produced from LI data of each of the three WISE epochs (see Table |Cl) shown in order from left to right, with the SN 
position marked by a green arrow. North is up and East left, and images are approximately 2'.5 on a side. 



galaxy results in non-uniform underlying residuals. The com- 
bination on an uneven background with the large WISE PSF 
may explain why the source was not picked up by the auto- 
mated pipeline. 

The Atlas image presented in Fig.[2]is made from a combi- 
nation of all images from the first epoch (2010 Jan) of WISE 
imaging and some images from the second epoch (2010 Aug). 
Comparison with our coadde d im ages produced from the Lib 
frames (see description in § C.2 and lower panels shown in 
Fig. [2}, we find that the signal in the Atlas images is com- 
pletely dominated by first epoch (2010 Jan). 

Subsequent observations by WISE (in 2010 Aug and 2011 
Jan) do not show any appreciable detections. Detection upper 
limits at epochs of ^300 and 470 days post-discovery in Ta- 
ble C2 show that the SN declined by a factor of at least 1.4 in 
flux (and probably more) in both bands on timescales of ^6 



months or less. 

3.5. Optical spectrum 

The Subaru spectrum of SN 2009js is shown in Fig. [4] 
The dereddened data have a blue spectral shape over the op- 
tical regime as expected from the hot temperature Tw7000 
K around 16 days post-discovery. A P Cygni profile of the 
absorbed blue wing for a strong emission line is most clearly 
visible for Ha, but also for H/3. Many other significant ab- 
sorption features can be discerned over the full wavelength 
range. 

The Na I D doublet near 5891 A is a tracer of line-of-sight 
extinction, and our spectrum shows two Na I absorption fea- 
tures: one in the restframe of the host galaxy, and one consis- 
tent with extinction arising in the Milky Way. The equivalent 
widths of the features are wl.l A and 3.6 A, respectively. The 
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FIG. 3. — Reddening-corrected BVRI light curves of SN 2009js compared to those of SN 2005cs (thin dotted curves) and SN 2008in (dashed) covering both 
the Kanata and NTT monitoring data to an epoch of about 360 days. Both Galactic and local reddening corrections have been applied for all objects. The cyan 
upper-limit arrows and detections on the left denote the pre-explosion limits and early time photometry from the literature. This is unfiltered photometry, closely 
matched to the V band (see Appendix for details). The circled symbols on the right are the NTT VRI magnitudes. Light curves for SN 2005cs and SN 2008in 
have been scaled in order to approximately match during plateau, for comparison. 



Milky Way feature is the stronger of the two, which is consis- 



tent with the reddening trend determined in j|3.1|and 3.2 



In order to classify the SN spectrum, we used the SNID code 
(Blondin & Tonry 2007 ) in order to compare this source to 
previous events . The returned clos est match is to the Type 
IIP SN 2005cs ( |Pastorello et al.|2009) . The optical spectrum 
of this source at a phase of ~17 days is overplotted in Fig. |4] 
and shows an excellent match to SN 2009js. Measuring the 
centroid wavelength of the trough of the Ha P Cygni absorp- 
tion component in the restframe spectrum yields a blueshift of 
5216±30 km s" 1 for SN 2009js. 

4. DISCUSSION 

4.1. Explosion epoch 

The exact explosion epoch of SN 2009js is unknown. Early- 
time discovery epoch photometry aligns well with our subse- 
quent Kanata monitoring meaning that we do not probe the 
initial flux rise. Using the non-detection at an epoch of -11 
days (Table |AT) i the explosion can be constrained to have oc- 
curred at MJD 55109.94 (2009 Oct 05.94) to within an uncer- 
tainty of 5.5 days. 

4.2. Evolution of temperature and bolometric luminosity 

Multi-band photometry is available for multiple epochs of 
Kanata monitoring, and for one late epoch at the NTT (Ta- 
ble [AT}. Each epoch was corrected for reddening due to the 
Galaxy and that local to the source, and then fitted with single 
Planck functions under the assumption that the early evolution 
is dominated by photospheric emission. Of the 25 nights of 



optical Kanata photometry, epochs +3.4 and +107.1 coincide 
closely with the MIR detections of SN 2009js with Spitzer 
and WISE, respectively. In these two cases, the optical and 
MIR data were fitted simultaneously. Best fit parameters and 
uncertainties were measured using the MPFIT package ( |Mark-| 
|wardt||2009| in idl. 15 This allows measurement of photo- 
spheric temperatures (T) and (quasi)bolometric luminosities 
(^ubvriX both of which are plotted as functions of time in 
Figs.|5]and|6]respectively. Three representative epoch fits are 
shown in Fig. [7] 

There is an apparent fast rise in T between the first two 
Kanata temperature measurements (a timespan of three days), 
which is intriguing given that none of the light curves in Fig. [3] 
show any dramatic increase at this phase. Some caution is 
necessary in interpretation of this rise, given that it is only a 
single measurement of temperature change and the uncertain- 
ties are somewhat large. The subsequent evolution shows a 
decline from a peak of ^8500 to <5000 K over a timespan of 
about 120 days. Data in all optical filters were not available 
for each epoch, resulting in large uncertainties when only a 
few filters were observed. But the systematic trend of temper- 
ature evolution is clear irrespective of these variations. 

Fig. [6] shows the change in Lubvri> the optical quasi- 
bolometric luminosity. Although we do not have U band 
monitoring, Ltjbvri can be measured as the integration under 
the fitted Planck curve between wavelengths of 0.35-1 /zm. 
This parametrization of the optical bolometric power has been 

15 http://www.exelisvis.com/ProductsServices/IDL.aspx 
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FIG. 4. — Subaru/FOCAS spectrum of SN 2009js at 1 6 days post-discovery (black) compared to the Mt. Ekar 1 .82m AFOSC spectrum of SN 2()05cs 17 days 
post- discovery (red cu rve observed on 2005-07-14; from Pastorello et al. 2009 1, and with SN 2008in 14 days post-explosion (blue curve observed on 2009-01-07; 
from |Roy et al.|201l) . The comparison spectra of SN 2005cs and SN 2008in were normalized to the continuum of SN 2009js around 7000 A, and then offset 
slightly for better legibility. The match between SN 2009js and SN 2005cs is excellent, except that the continuum slope of SN 2005cs is slightly bluer as a result 
of its hotter temperature. SN 2008in was further scaled in order to match the peak strength of the Ha emission line. SN 2008in shows a significant difference 
with respect to the other two objects in terms of having a much shallower blueshifted Ha absorption feature. 



chosen to facilitate direct comparisons with other measure- 
ments in the literature, which will be discussed shortly. The 
peak luminosity is log Ltjbvri ( er g s" 1 )^!^. This quasi- 
bolometric light curve also shows four distinct phases: 

1. an initial fast fade within the first ~20 days; 

2. a period of approximately constant output for the next 
^70 days; 

3. then a sharp decline in the transition between the 
plateau and nebular phases to the last probed Kanata 
datapoint; 

4. and finally, an apparently gentler decline between the 
Kanata and NTT epochs. 

The mean change in -Lubvri during these four periods are 
approximately: 1) -0.01 1; 2) -0.001; 3) -0.024 and 4) -0.006, 
as measured in log luminosity (dex) per day. 

The total bolometric power (Lboi) may be computed by 
integrating under the full Planck curve, including contribu- 
tions from the far-ultraviolet and infrared regimes. This is a 
factor of 1.50 higher than £ubvri> on average. Integrating 
£bo1 over all epochs returns a total radiative energy of about 
5.2xl0 48 erg. 

4.3. SN 2009js: comparisons with SN 2005cs and SN 2008in 

The optical light curves discussed above (especially in the 
R and / bands, as well as in -Ltjbvri) show the presence of a 
clear plateau phase. This confirms the Type IIP nature of the 
source, initial suggestions of which were based upon very- 
early-time spectroscopy (Silverman et al. 2009). The long- 
term VRI light curves shown in Fig.[3]suggest that the change 
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FIG. 5. — Evolution in blackbody temperature after correction for local and 
Galactic reddening. 



in brightness from mid-plateau to late nebular phase is overall 
similar to that of the subluminous SN 2005cs. 

Silver man 



Based up on early-time optical spectroscopy, 
et al. ( 2009 [ I measured an Ha absorption redshift of about 
7200 km s _i at two days after maximum brightness, and sug- 
gested a match to SN 2005cs. Peak brightness corresponds to 
a phase of 2 days according to |Pastorello et al.] ([2009). Our 
measurement of the blueshift in the Subaru spectrum obtained 
at a pha se o f about 16 days post-discovery is 2000 km s _1 
lower (§ 3.5 1, as expected from the gradual lowering with time 
of the outer photospheric layer that is optically-thick. This 
Subaru spectrum also displays remarkable affinity to that of 
SN 2005cs at 17 days (Fig. HI). Although there is some uncer- 
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FIG. 6. — Evolution in quasi-bolometric power (integrated under the Planck function between 3500A and 1 fim) after correction for local and Galactic 
reddening. Data for comparison SNe are taken from Fig. 9 of Pastorello et al. 1 2009 1 except for SN 2008in, for which Ltjbvri i s computed from the photometric 
data presented by |Roy et al. 1(2011) in a manner identical to that tor SN 20U9js. 
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0.1 1 10 0.1 1 10 0.1 1 10 

jim |im |im 



FIG. 7. — Blackbody fits to multi-wavelength epoch SEDs, corrected for local and Galactic reddening. Three epochs close to the beginning (epoch +3.4), 
middle (epoch +56.2) and end (epoch +107.1) of Kanata monitoring are shown for illustration. The first and last of these include MIR detections. 



tainty in determination of the explosion epoch of SN 2009js, 
these matches provide evidence (in addition to that discussed 



4.1 1 that SN 2009js was discovered close in time of ex- 



plosion (i.e. within a few days). All these lines of evidence 
suggest that SN 2009js was a 2005cs-like event. 

Both sources show peak optical bolometric luminosities 
of well below 10 42 erg s _1 , above which 'normal' Type II 
SNe generally lie. Furthermore, both sources have absolute 
F-band magnitudes much lower (intrinsically fainter) than 
My=-16. This appears to support a subluminous classifica- 
tion for SN 2009js also. But Fig. [6] shows that the plateau 
phase luminosity of SN 2009js is instead much more closely 
matched to that of another source, SN 2008in. This w as also 
a Type IIP event, studied in detail by |Roy et al. ( 201 1| >. It has 
been classified as an intermediate luminosity source, bridg- 
ing the 'gap' between the subluminous and normal popula- 
tions. We computed LuBVRi(2008in) by using the multi-band 



photometry and extinction determination presented by Roy 



|et al.| ( |201 1) , by using Planck-function fitting exactly as for 
SN 2009js. Almost identical results are obtained if we use the 
quasi-bolometric lightcurves (in terms of UVOIR) presented 
in Fig. 12 of Roy et al. ( 20 1 1 ] > and use a mean conversion 



factor to the U BVRI wavelength range over the lightcurve. 
In addition, Fig. 8] shows single-band absolute magnitude 
light curves, avoiding any uncertainties related to bolometric 
corrections, and again suggests a reasonable match between 
SN 2008in and SN 2009js during the plateau phase. On the 
other hand, the main difference with respect to SN 2008in is 
in terms of the spectra: Fig.|4]shows that significantly broader 
Ha in emission for SN 2008in and a blueshifted absorption 
feature which is weaker by a factor of about 2 (in terms of 
equivalent w idth) as compared to the other two sources. Roy 
et al. ( 201 1) have also studied the evolution of the spectral ve- 



locities for SN 2008in and found significantly larger velocities 
than for SN 2005cs over the first 100 days. 

Thus, it appears that SN 2009js shares mixed characteristics 
of both 'subluminous' and 'intermediate' luminosity events. 
We discuss other similarities and differences between these 
three sources in the following sections. Ultimately, however, 
the exact nomenclature of the assigned class is not crucial. 
Rather, the main question is whether the underlying progen- 
itor and explosion properties of these events are significantly 
different or not. 
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4.3. 1 . Further detailed comparisons 

1. The plateau luminosity of SN 2009js is higher than 
that of SN 2005cs by a factor of about 2, and more closely 
matched to SN 2008in, as mentioned in the previous sec- 
tion. During transition between plateau and nebular phase, 
however, SN 2009js fades and falls below the luminosity of 
SN 2005cs. During the late nebular phase, the luminosity of 
SN 2009js lies in between that of SN 2005cs and SN 2008in. 
This is readily apparent from the quasi-bolometric luminosity 
light curve and the absolute V band light curve in Figs.[6]and 

m 

2. The plateau phase is somewhat shorter than that of 
SN 2005cs, which is why the luminosity of SN 2009js drops 
below that of SN 2005cs during the transition phase. In Eq. 
2, we have measured ipT(2009js)=l 1 1 days. We carried out 
an identical fit for the light curve of SN 2005cs and found a 
significantly larger £pt(2005cs)=130.8±0.2 days. The differ- 
ence of about 20 d ays is larger than the uncertainties on the 
explosion epoch (§ 4.1 1. A comparison with other SNe yields 
a similar result, a fact that is best visible in Fig. [9] This figure 
shows the evolution of Lubvri values for various SNe as in 
Fig. [6] with the only change being that all sources are shifted 
vertically (normalized) in order to approximately match the 
early-time peak of SN 2009js. Apparently, the drop in L-q \ 
for SN 2009js occurs from a few days before, to up to two 
weeks before, that in most other sources: SNe 2005cs, 2004et, 
1999em and 1999gi. However, the plateau length does appear 
to match that of SN 2008in, which has also been measured to 
be short 
inFi g.[3]l. 



days according to Roy et al. 201 1 also apparent 



3. |Pastorello et al. ( 2006[ ) found unusually blue U-B colors 
for SN 2005cs, corresponding to very high continuum tem- 
peratures T of up to 3 x 10 4 K at early times, which cooled to 
T~7500 K at 35 days. According to our multi-epoch black- 
body fits (§ 4.2 1, T(2009js) is significantly lower over the en- 
tire period probed (Fig.[6]i. Fig.[5]shows a reasonable match to 
the temperature values found for SN 1999br and SN 2008in. 
We do not have U band monitoring data for SN 2009js, so one 
may question whether we miss a significant fraction of the 
flux and underestimate the photospheric temperature. How- 
ever, our multi-band optical photometry appears to satisfac- 
torily describe the Wien tail of the optical data (see, e.g., the 
decrease towards shorter wavelengths in Fig. |7J. So we are 
unlikely to be missing any significant component in the U 
band, and T(2009js) is then significantly lower than that of 
SN 2005cs. This conclusion is robust to uncertai nties i n hos t 
galaxy reddening, which is estimated to be low (§ 3.2 § |3.5[ ). 

4. Fig. [9] shows that SN 2005cs underwent a larger lumi- 
nosity decrement between peak and plateau, as compared to 
SN 2009js. In order to quantify this, we define a new parame- 
ter Alog£ as the difference between the peak luminosity and 
that measured at 50 days post-peak in log luminosity units: 



Alog£ = log£ (tpoak) - logL (Vak+50d) 



(7) 



We have Alog£(2009js)=0.16 dex and 
Alog£(2005cs)=0.36 dex, when measuring the lumi- 
nosity in terms of Lubvri- The match is even worse when 
compared to SN 2008in, for which Alog£(2008in)=0.47 dex. 
This indicates a comparatively lower amount of radiative 
cooling for SN 2009js during adiabatic expansion following 
shock breakout. In this respect, SN 2009js is closer matched 
to SN 1999em (Alog£=0.19) and SN 1999gi (Alog£=0.22). 



4.4. Nebular phase 

The NTT data reveal the late phase evolution of SN 2009js. 
The VRI light curves in terms of observed flux (Fig. [3]) show 
that the relative decline from mid-plateau to the nebular phase 
is very similar to that of SN 2005cs in all three available fil- 
ters. Although we do not probe the beginning of the nebular 
tail in our Kanata monitoring, the comparison with SN 2005cs 
implies that the I band is likely to have covered the bulk of 
transition phase after the plateau. This allowed us to mea- 
sure the slope of flux decline in the nebular phase in § 3.3 as 
Po=0.009 mag day -1 , a value which lies within the range of 



values found for other sources in the compilation of Olivares 
[eTaLl ( |20T0| . 

The nebular phase is the result of energy release dur- 
ing radioactive decay of 56 Ni. Thus, the luminosity in 
this phase is a good tracer of the amount of 56 Ni pro- 
duced. The fitted blackbody temperature in this phase is 
3530±1000 K. In terms of (quasi)bolometric optical power, 
£uBVRi(2009js)=2.08(±0.59)xl0 39 erg s" 1 . The luminos- 
ity uncertainty was determined by Monte Carlo sampling of 
fluxes drawn from a normal distribution with standard de- 
viation equal to the EFOSC photometric errors in each fil- 
ter, followed by integration of the fitted Planck curves to 
give an ensemble of randomized luminosities. This lumi- 
nosity is 2.4(±0.7)xi UB vRi(2005cs) at an epoch of +360 
days, and about 0.4(±0.1)x_LuBVRi(2008in); in other words, 
SN 2009js lies between the other tw o SNe. The 56 Ni mass 
produced in SN 2008in according to |Roy et aL] ( |201 1[ ) was 
0.015 M . Scaling this by luminosity implies MNi(2009js)« 

0. 006(±0.002) M Q . Extrapolating instead from the 56 Ni 
mass synthesized in SN 2005cs (~0.003 M . . according 
to |Pastorello et al. 2009) > produces a very close result, 

1. e. iW Ni (2009js)«0.007(±0.002) M . Assuming a 50% 
higher bolometric correction (see § 4.2 1, this may increase to 
M Ni (2009js) wO.01 1(±0.003) M . 

An important caveat regarding measurement of the lumi- 
nosity is that a Planck function may be a poor approxima- 
tion to the source spectrum during the nebular phase. This 
has been investigated by Hamuy (2001) who suggests that al- 
though such an approximation is worse in the nebular phase, 
the Planck function still gives comparable results to direct 
integration over the full optical-infrared regime, to within 
15% accuracy. Using the bolometric correction of 0.26 mag 
from the ^-band found by |Hamuy| pOOT} for SN 1987 A 
and SN 1999em, we compute a nebular phase total bolomet- 
ric luminosity (LboO of 3.6(±2.7)x 10 39 erg s -1 which is 
only 1.7x larger than (and consistent within the reddening- 
corrected EFOSC F-band flux uncertainty with) the value 
of Lubvri use d above and plotted in Fig. [6] If we instead 
use the R band absolute magnitude in order to avoid un- 
certainties related to bolometric corrections, then we have 
M R (2009js)wM R (2005cs)-0.20 (with a 40% reddening- 
corrected flux uncertainty) at an epoch of +360 days. Scaling 
M Ni then implies M Ni (2009js)«0.004(±0.0016) M . 

In summary, the above range of Mni is much smaller than 
synthesized in normal SNe II (for which Mni is generally 
greater than 0.05 M ), supporting the association with the 
subluminous class of events. 

4.5. Infrared detections: photosphere or dust emission? 

The Spitzer detections of SN 2009js at an epoch of 2 days 
post-discovery, or 7.5±5.5 days post-explosion (§ 3.4 1 repre- 



sent some of the earliest mid-infrared (MIR) SN detections, 
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FIG. 9. — Evolution in quasi-bolometric power (integrated under the Planck 
function between 3500A and 1 /im) after correction for local and Galactic 
reddening. 



especially beyond 4/im. Only SN 1987 A has comparably- 
early mid-infrared follow -up in the shorter-wavelength L 
band ( |Menzies et al.|1987) . Subluminous SNe have not been 
studied in detail in the mid-infrared so far. 16 The fact that mid- 
infrared detections of SN 2009js with both Spitzer and WISE 
were serendipitous (i.e. not part of any SN survey) empha- 
sizes the importance of examining archival data and all-sky 
surveys in detail. 

The question of dust production in SNe is a very active topic 
of research. Dust has been detected in several SNe using mid- 
and far-infrared observations ( e.g. [W ooden et al. 1993 Kotak 
et al.|2 009HMerkle et al |201 l[|Matsuura et al.|201 lt|Fox et aL 
Szalai & Vinko 2012). The pace of dust production 



20111 



16 Archival Spitzer data of SN 2005cs was studied by jSzalai & Vinko| 
(20121, but no reliable detection was found. 



in core-collapse SNe is uncertain, but is generally regarded 
to have a characteri stic timescale of seve ral months to one 
year or more (e.g. Nozawa et al. 2003 i. Early-time MIR 



excesses have been detected in some cases (e.g. in SN 1987 A 
by |Wooden et al.|[T993] and in SN 2004dj by |Meikle et al] 



20 



1). In SN 1987A, the excess is attributed to heating of 



circumstellar dust, while in SN 2004dj, it is instead associated 
with echos from newly formed dust within a cool dense shell. 

With only two MIR bands and the absence of any near- 
infrared monitoring, it is difficult to assess the contribution of 
dust to the spectral energy distributions (SEDs) of SN 2009js. 
However, we note that the spectral slopes between the two 
available mid-infrared bands are flatter than those expected 
from a blackbody fitted to the broadband (optical+MIR) data. 
This is easily visible in Fig.|7]and is found to be true for both 
Kanata epochs +3.4 days and +107.1 days when MIR detec- 
tions are available. In the former, the slope is only marginally 
flatter than the best-fitting Planck function within the photo- 
metric uncertainties. In the latter, on the other hand, the 4.6 
WISE W2 excess appears to be significantly larger. This 
may be consistent with the expectation of increased dust pro- 
duction on timescales of months. 

Although accurate assessment of the dust contribution is 
difficult, it is possible to compare SN 2009js with another 
Type IIP SN 2004dj, in which MIR e xcesses attributed to 
freshly formed dust have been observed ( |Meikle et al.|201 1\ . 



Fig. 10 compares the optical and MIR SED of these two 
sources at a similar epoch of 106 days. For the latter source, 
the data and triple blackbody fit are both taken from Meikle 
et al. ( |201 l| l and are normalized to match SN 2009js in the 



optical. From the fact that the MIR fluxes of SN 2009js lie 
well below those of SN 2004dj (normalized to the optical 
photospheric emission), it is immediately clear that dust heat- 
ing in SN 2009js is much less efficient than in SN 2004dj. 
In order to construct a comparative model of any potential 
dust emission in SN 2009js, we assumed that the photo- 
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spheres of the two sources appear identical at this epoch, i.e. 
both h ave a temperature T p i lo tosphere=7000 K (Meikle et al. 
201 1\ . Two blackbody functions were used to represent the 
photosphere and the dust respectively. The fitted model is 
shown in Fig. [10| and yields T dust (2009js)«630 K, which 
is lower by a factor of about 3 than T dust (2004dj)=1800 K 



from Meikle et al. ( 201 l| l. Assuming that the emission re- 
gion is optically th in and that grain c omposition is mainly 
amorphous carbon (Zubko et al. 1996| l, we find a dust mass 
of M dust (2OO9js)~3.OxlO-°M0. 1 ' 

The integrated dust luminosity is only about 3 % of the in- 
tegrated photospheric luminosity in SN 200 9js, whereas it i s 
about 15 % for SN 2004dj (see Table 3 of |Meikle et al.|201 1) . 
The absolute luminosity of the dust component for SN 2004dj 
is also higher, at about 2.5 times the value that we find for 
SN 2009js, even though the luminosity of the photosphere 
of SN 2004dj is abo ut half that of SN 2009js at this epoch. 

attributed the MIR excess to freshly- 
If the same were to 



Meikl e et al] pOTT 
formed dust within a cool dense shell 



be true in the case of SN 2009js and if both events had sim- 
ilar progenitors, this would imply much less dust production 
in SN 2009js. But the present data cannot distinguish freshly- 
formed dust from a scenario where the infrared excess origi- 
nates as an echo from pre-existing dust. This latter possibility 
may be more likely, as various statistical studies of SNe dust 
have indicated ( |Fox et al.|20TT]|Szalai & Vink6|20T2"l ). MIR 
detections at epochs of about 100 days may be too early for 
the infrared excess to be interpreted as thermal radiation from 
dust formed in the ejecta (in Type IIP SNe, dust gener ally 
condenses after about 300 days; e.g., Nozawa et al.|2003 1 



One other possibility to consider is that the 3.4 fj,m emis- 
sion is photospheric and not dust-related (as suggested in the 
final panel of Fig. [7]), whereas the 4.6 /im excess is the result 
of fundamental CO molecular line emission. Such molecular 
emission is thought to be an important precursor signal to the 
imminent onset of dust condensation in SNe eje cta, and has 
been discussed in the context of several SNe (e.g., Kotak et al. 
[20061 |Meikle ft al.|2011) . 

To summarize, the comparison with SN 2004dj suggests 
evidence of heated dust emission in SN 2009js at an epoch 
of about 107 days. Alternatively, the MIR excess above the 
photosphere may be a sign of dust-precursor molecular emis- 
sion. Both of these scenarios are interesting, and it is impor- 
tant to extend studies of dust formation to the subluminous 
SNe regime with further detailed follow-up work. 

4.6. The nature of the progenitor star and the SN explosion 

The nature of the progenitor of SN 2009js is presently un- 
known. The fact that the source is 2005cs-like may suggest 
a similar progenitor, which is constrained to be a moderate- 
mass (8-15 M Q ) star (jPastorello et al. 2009 ). Progenitors of 
subluminous events may generally lie at the lowe r mass range 
of stars that produce core-collapse events (e.g. IFraser et al. 



201 1 1. The mass constraint on the intermediate luminosity 
source SN 2008in is obtained only from analytic modeling of 
the light cu rve and spectra: an upper-limit of 20 M Q (Roy 
|et al.|2011) . 

We have Spitzer pre-explosion images for SN 2009js, but 

17 A predominantly Silicate dust composition, or a much cooler dust com- 
ponent, cannot be ruled out, but these are less likely because a prominent 
detection in the WISE W3 filter (encompassing the Silicate emission feature 
wavelength) is predicted in these cases at a level of ~ 1 mJy. This flux level 
corresponds closely to the canonical flux limit of the WISE all-sky survey 
I Wright et al. 20101, and there is no such detection in the all-sky release. 




Wavelength (|j.m) 

FIG. 10.— Comparison of the optical and MIR SEDs of SN 2009js with 
SN 2004dj. Two-blackbody fit to the optical and mid-IR data of SN 2009js 
is denoted by the black continuous curves and filled circles at an epoch of 
~107 days (Table [AT) . The photospheric temperatures of the two model 
sources have been matched. For a Carbon dust model, the mass of the dust 
fitting the mid-IR data is ~3x 10 Mq. In light color and unfilled green 
squares, the SE P and triple blackbo dy fit for SN 2004dj at an epoch of 106 
days are shown (Meikle et al. 2011 1. These have been normalized to match 
SN 2009js in the optical. 

these do not provide stringent constraints on the progenitor. 
We instead use established analytical relations in order to ap- 
proximately estimate the explosion energy (E), ejected mass 
(M e j) and the pre-S N stellar radius (i? pre _ sN)- For this pur- 
pose, Eqs^l -3 of |Litvinova & Nadezhin| ( |1985| see also 
Popov 1993) may be used, with the caveats that these rela- 



tions ignore some important effects such as nickel heating 
and line opacities. The required observables are the plateau 
length (At) and the absolute V magnitude Aly at mid-plateau, 
and the photospheric velocity (u pfl ) also at mid-plateau. For 
SN 2009js, we have Ai«lll±5.5 days. This is taken to be 
tpT (Eq. 2), assuming that the plateau begins after the ini- 
tial rise of a few days (which we do not probe), and also in- 
cludes the uncertainty on the determination of the explosion 
epoch. 7\/v=-15.2±0.3 at an epoch of 55 days. The absence 
of spectroscopic monitoring means that « pn at mid-plateau is 
unkn own . But the good spectral template match to SN 2005cs 
(cf. § 4.3 i suggests that using an identical velocity evolution is 



a reasonable way to proceed. Interpolating the Sc II velocity 
evolution for SN 2005cs found by |Pastorello et al.| ( [2009] > at an 
epoch of about +55 days, we find u p h~1600 km s to which 
we assign a generous systematic uncertainty of 400 km s _1 to 
account for the absence of observations for SN 2009js. 18 This 
yields 



E = 0.14 ±0.11 foe 



3.9 ± 4.8 M, 







i? pre _ SN = 415±340 R C: 



(8) 
(9) 
(10) 



The ejecta mass is about 9 M . Adding about 1 M Q for 
stellar wind mass losses and another 1 .4 M for the compact 
remnant assumed to collapse during the SN, the initial pro- 
genitor mas s is estimated to be about 11 ±5 M Q . |Maguire| 
et al. (2010) have shown that progenitor mass estimates based 
upon analytical relations are generally higher than those based 

18 Such an uncert ainty would also encompass the Fe II velocity observed 
by |Roy et al.| {20l 1 1 for SN 2008in at a similar epoch, at the 2a confidence 
level! 
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upon pre-explosion images or detection limits, and some addi- 
tional caution is thus required when using these values. How- 
ever, a useful statement which can be made is that these es- 
timates of M e j and Mp rogcn it r are consistent with those de- 



rived for_SN_2005csusing various methods (Past orello et al.| 
2009{ paguire et al.|2010") . 

The explosion energy derived above for SN 2009js is lower 
than that for SN 2008in (0.5 foe; |Roy et al.|[20TT) . The 
val ue of E derived simila rly for SN 2005cs (0.17±0.08 foe) 
by Maguire et al. (2010), on the other hand, is consistent 
with that of SN 2009js, using identical analytical estimates 
as above. These results provide additional support of a close 
association between SN 2009js and SN 2005cs. 19 

5. SUMMARY 

Early-time photometry previously reported constrains the 
explosion epoch of SN 2009js to MJD 55109.94±5.5 days 

(§ED- 

Kanata optical monitoring shows S N 20 09js to have charac- 
teristics typical for a Type IIP event (§ 4.2 1. The characteristic 



plateau length (as measured by the mid-point of th e tra nsition 
between plateau and nebular phases) is 1 1 1 days (§ 3.3 1. Line- 
of-sight extinction is found to be Ay(Galaxy)gaQ.95 mags 



(§ 3.1 1 and A v (host galaxy)=0.18±0.38 mags (§ 3.2 1, i.e. the 



reddening local to the SN is lower than that due to our Galaxy. 

Subaru optical spectroscopy shows a very good match 
to the spectrum of the well-studied subluminous Type IIP 
SN 2005cs (§ [33J. These facts suggest that SN 2009js was 



a SN 2005cs-like event (§|43i. 

We are able to chart the evolution of the photospheric tem- 
perature and optical bolometric power of SN 2009js (§ 4.2 1 



The plateau-phase luminosity of SN 2009js is higher than that 
of SN 2005cs, but significantly lower than normal Type IIP 
events such as SN 1999em (Fig. [6]) and SN 1999gi (Fig. [8}. 
Some ot her di fferences with respect to SN 2005cs are also 
noted (§ 



4.3.1 



The initial post-peak adiabatic phase cool- 
ing is less prominent than in SN 2005cs, as quantified by 
a new parameter Alog£ measuring the difference in log lu- 
minosities at peak, and 50 days post-peak. Moreover, the 
plateau length of SN 2009js is significantly shorter than in 
SN 2005cs, and more similar to that of the intermediate lumi- 
nosity event SN 2008in. The absolute mid-plateau luminosity 
of SN 2008in also closely matches that of SN 2009js. On the 
other hand, SN 2008in shows a significantly higher nebular 
luminosity, as well as a larger cooling dip Alog£. In addi- 
tion, its photospheric velocity determined from spectroscopy 
is significantly higher than SN 2005cs and SN 2009js. These 
facts imply that SN 2009js has characteristics common to both 
subluminous and intermediate luminosity events. 

The mass of synt hesi zed radioactive 56 Ni is found to be 
0.004-0.011 M Q (§|4~4]l. This is similar to the synthesis in 



Using known analytical relations, a low explosion en- 
ergy _E=0.14±0.11 foe is derived. This is lower than for 
SN 2008in, within the caveats of use of the analytical rela- 
tions. The constraints on the ejecta mass (8.9±4.8 M Q ) and 
progenitor mass, as well as on the explosion energy, are fully 
consistent with those derived for SN 2005cs using similar 
methods (§ |4.6| l. 

SN 2009js is the first subluminous (or intermediate lumi- 
nosity) SN to be studied in the mid-infrared (MIR). Spitzer 
detected SN 2009js serendipitously at a very early epoch, and 
WISE monitored the field on three occasions to an epoch of 
about +470 days. A significant MIR excess above the pho- 
tosphere is found at 4.6 pm at an epoch of ^107 days post- 
discovery. A tentative dust model fit implies a mass of amor- 
phous Carbonaceous dust of about 3xl0 -5 M Q , and a low 
efficiency of dust heating as compared to SN 2004dj (§ |4.5| ). 
Alternatively, this may be a sign of fundamental CO molecu- 
lar emission which is a precursor to dust formation. 

Low luminosity events constitute only a few per cent of the 
Type II SN population (e.g. Pastorello et al. 2004), and few 
studies exist on how these connect to the 'normal' population. 
Our observations of SN 2009js thus add to the list of sources 
in this class with relatively well-sampled multi-wavelength 
light curves and provide first insight into their mid-infrared 
properties. 
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keeping the caveat in mind that the E estimate is quite sensitive to u p ^, 
for which we do not have a direct measurement. 



APPENDIX 
OPTICAL PHOTOMETRY 

Kanata optical imaging 

The 1.5 m Kanata telescope ( |Uemura & Kan ata Team 2009), part of the Higashi-Hiroshima Observatory at an altitude of 503-m 
above sea level, is used by our group for regular monitoring of SNe. HOWPol BVRqIq imaging provides a 15'-diameter field 
of view at a pixel scale of 0'.'29. 

Monitoring started on 2009 Oct 14.8 at an epoch of +3.4 days and continued for about 4 months. A total of 25 epochs of data 
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on different nights were obtained during this time. 

The images were reduced according to normal procedures for CCD photometry. We performed point-spread-function fitting 
photometry using the DAOPHOT package in IRAF 20 . We also observed PG 0231+051 as a standard star field from the Landolt 
standard star catalog (Landolt 1992 ) when the night was photometric. The magnitudes of local comparison stars 21 and the color 
terms were calibrated using these standard stars. Photometric uncertainties include Poisson measurement error and standard 
deviation resulting from calibration against different comparison stars. Zeropoints are those from the canonical filter definitions 
of ( [Bessell et al.|1998) . 



NTT late-time optical imaging 

NTT/EFOSC was used for observing SN 2009js on 2010 Oct 6. 22 Observing conditions were cloudy with a median seeing of 
about 1'.'9 in the V band. A series of five consecutive exposures (100 s each, giving a total exposure time per filter of 500 s) were 
taken in VRI starting at UT05:53 at a median airmass of 1.48. 

Reduction was carried out using the EFOSC2 pipeline and standard procedures for analysis ( |Izzo|201 T} . The reduced expo- 
sures were combined to a cleaned final image. Since the source is relatively faint at this late epoch, the background was carefully 
measured using adjacent apertures that include emission from the sky as well as the underlying host galaxy. No photometric stan- 
dard star observations were obtained that night. We therefore carried out relative photometry with respect to nearby comparison 
stars, and then tied the absolute calibration to the Kanata data using stars common to the Kanata and NTT images. This also 
alleviates the issue of small differences in filter response between the VRI filters of the two telescopes. The resultant photometry 
is listed in Table I All 



Early filterless photometry 

In addition to the abov e datasets analyzed by us, we also made direct use of pre-discovery and discovery photometry reported 
in CBET 1969 (Nakano et al. 2009; Silverman et al.|2009). These data were obtained in filterless CCD mode and are listed at the 



top of Table Al Although unfiltered CCDs have broad response functions, their effective central wavelengths are typically close 
to that of the V band. |Riess et aT (|T999 ) have determined cross-calibration equations between unfiltered and V band photometry, 
based upon the B — V color. Although we do not know the SN colors at the discovery epoch, the first Kanata measurements 
corrected for reddening (see §|3.1|and|3.2|i imply B — F«0.3, as expected for a hot blackbody characterizing the early-time SN 
spectrum. Using Eq. 3 of Riess et al. ( 1999 ) and the transformation constant from Table 2 therein implies that the unfiltered mag 
is expected to be within 0. 1 mags of the V-band mag. We thus use the unfiltered mags directly as a proxy for the V band, with a 
systematic uncertainty of 10%. 



OPTICAL SPECTROSCOPY 

SN 2009js was observed on UT 2009-10-27 (MJD = 5 5131 .6), about 16 days p ost-discovery with the Subaru telescope equipped 
with Faint Object Camera and Spectrograph (FOCAS; [Kashikawa et aT1|2002| >. The blue (3800-7000 A) and red (6000-10000 
A) parts of the spectrum were taken separately. The blue part was taken with the 300B grism without an order-cut filter, while 
the red part was taken with the 300R grism with the 058 filter. These configurations give a wavelength resolution of 10 A (500 
km s _1 near the center of the spectrum, R = A/AA ~ 600). For both settings, we used a center slit of 0V 8 width. The target was 
integrated for 600 s in both settings. The typical seeing during the observations was O'.'S in the V band, with median airmass of 
about 1.65. Reduction was carried out in IRAF. Wavelength calibration was performed using the Th-Ar arc lamp. Flux calibration 
was carried out using spectrophotometric star observations on the same night. Finally, the blue and red part of the spectrum were 
combined. 



INFRARED PHOTOMETRY 
Spitzer mid-infrared imaging 

Archival MIR imaging data of NGC 918 are available from the Spitzer archive in both bands IRAC1 and IRAC2. These were 
taken in late 2009, a few months into the warm mission phase. The instrument point response functions are complex, but imaging 
quality is better than 2" in terms of equivalent Gaussian full- width at half-maximum measurements 23 . Level 2 post-BCD IRAC 
mapping mosaics from the post-cryogenic phase of the mission, processed by standard pipeline software version S 18. 12.0, were 
downloaded and directly used for the present analysis. 

Spitzer happened to observe the galaxy at epochs of about -35 days and +2 days. For each of the IRAC1 and IRAC2 bands, 
the pre-explosion epoch image was subtracted from the post-discovery one to remove background structured emission from the 
host galaxy. Source photometry was carried out in this subtracted image using a large aperture to encompass the total flux. 
Conversion from image units to final fluxes was done using standard procedures and conversion factors 24 . As a cross-check, we 



also carried out detection of all sources within a few arcmin of the host galaxy nucleus using the SEXtractor package ( Bertin & 
Arnouts|1996 i and computed a zeropoint for each image based upon cross-calibration with objects common to the WISE all-sky 



catalog for which all magnitudes are publicly-available. The bandpass difference between IRAC1 and WISE Wl is small and was 

20 http://iraf.noao.edu 

21 These stars have the following IDs in the USNO-B1.0 survey jZacharias et al |2004| : 1084-0033855, 1084-0033919 and 1085-0027773. 

22 Based on observations obtained during ESO program 1 84. D-l 140(E), P.l. S. Benetti. 

23 http://irsa.ipac.caltech.edU/data/SPITZER/docs/irac/iracinstrumenthandbook/5 

24 http://irsa.ipac.caltech.edu/data/SPITZER/docs/irac/iracinstrumenthandbook/IRACJnstrument_Handbook.pdf 
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TABLE A 1 
Optical photometry 
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EFOSC 

23.45±0.63 (0.004±0.003) 21.95±0.30 (0.010±0.004) 21.10±0.10 (0.014±0.003) 



^Denotes days post-discovery (2009 Oct 11.44). 
For the unfiltered early-time photometry, the observers are Itagaki (denoted by superscript /), KAIT (superscript K) and Yusa (superscript Y). The 
measurement by Yusa is not included in our analysis because it is brighter by 0.5 mags as compared to all other photometry obtained both before and after. No 

uncertainties are reported and it is impossible to judge its veracity. 
Listed magnitudes and magnitude limits are the observed ones in the B, V, R and / bands, whereas fluxes have been corrected for Galactic and host reddening. 

Errors in flux (lcr) include corresponding systematic dereddening uncertainties. 



ignored ( Jarrett et al.|2011 1. Similarly, IRAC2 and WISE W2 bandpass differences were ignored. We found our cross-check to 
be consistent with the direct Spitzer photometry. 

The log of observations is listed in Table CI images are presented in Fig. [T]and photometry tabulated in Table C2 Non- 
detection upper-limits in the pre-explosion images were determined by the addition of artificial point sources using the addstar 
routine in IRAF. This procedure involved generation of a point spread function (PSF) from nearby stars of known magnitude, 
followed by an iterative procedure of injecting scaled versions of the PSF into the coadded images at the position of the SN and 
running automatic source detection using the S EXTRACTOR package ( |Bertin & Arnouts|1996| l to determine the faintest reliable 
artificial source detection for which we assumed a 2er detection threshold. 

WISE mid-infrared imaging 

The Wide-field Infrared Survey Explorer (WISE) satellite (Wri ght et al.|2 010) carried out sky surveys in four bands (W1-W4 
centered on wavelengths of «3.4, 4.6, 12 and 22 fim, respectively), at an effective angular resolution of ss6" in W1-W3, and 12" 
in W4. Based upon the cryogen status, the publicly-released data products are divided into 1) an all-sky release, 2) a 3-band cryo 
release, and 3) a NEOWISE post-cryo release. Combined, these cover the period of 2010 Jan 7 to 2011 Feb l 25 . The sensitivity 
in the two short wavelength bands (Wl and W2) did not change dramatically with time, whereas the nominal sensitivities in W3 



' http://wise2.ipac.caltech.edu/docs/release/allsky/expsup/ 
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and W4 is degraded in the latter two releases. As of mid-2012, Atlas images as well as catalogs of measured positions, fluxes and 
other parameters of pipeline-detected sources are available for the all-sky release. These are the product of coaddition of several 
frames and have undergone multiple levels of calibration; they are referred to as Level 3 (L3) products. For the post-cryo data, 
only a preliminary release is available, including individual single -pass calibrated (level 1; Lib) data. SN 2009js is not included 
in the catalog of pipeline-detected sources of any release. 26 
Available WISE images for the field of NGC 918 cover three main epochs, starting from about 107 days after the discovery of 



SN 2009js. Subsequent epochs are separated by approximately six months. These are detailed in Table CI No pre-explosion 



WISE image for subtraction is available, so photometry at the position of the SN was carried out directly on Atlas and post-cryo 
images covering the SN position. 

For the fully-calibrated Atlas products, we used procedures recommended by the WISE team. 27 Small aperture photometry was 
carried out at the SN position as well as on bright nearby field stars with cataloged magnitudes used for relative calibration. The 
background was carefully selected using multiple adjacent apertures in order to avoid other point sources while encompassing 
local host galaxy emission. Aperture corrections, as well as correlated and confusion noise estimation for uncertainty measure- 
ments were included according to the instructions in the WISE explanatory supplement 28 . Magnitudes were converted to fluxes 
using standard zeropoints (Wright et a l . |20 1 0] > . 

The publicly-released Atlas (L3) image is made by coaddition of multiple exposures obtained in 2010 Jan, plus some exposures 
taken in 2010 Aug. We also produced our own coadded images for each of the three main epochs, starting from the Lib frames. 
First, the Lib data were sieved by rem oval of low quality frames affected by elevated background due to Moon proximity and 
other artifacts 29 . The swarp package ( jBertin et al.||2002 i was then used for background subtracting, resampling, shifting and 
median combining the good frames. Upper limits were determined identically as with Spitzer. 

Final images from the various WISE epochs are shown in Fig.|2]and photometry is listed in Table C2 



TABLE Cl 
Infrared Observation log 



Date (MJD) 
UT 


Telescope/Instrument/Survey 


Obs. 


identifier 


Bands 


Exposure time 

s 


2009 Sep 06.7 (55080.7) 

2009 Oct 13.6(55117.6) 

2010 Jan 26.1 (55222. l)t 

2010 Aug 02.8 (55410.8)+ 

2011 Jan 25.7 (55586.7)t 


Spitzer/IRAC/SAG 
Spfe_r/IRAC/S4G 
WISE/ All-sky release 
W/S£/All-sky release 
WYSS/Post-cryo release 


AOR= 
AOR= 


=30590976 
=30606848 


3.6 A*m (JRAC1), 4.5 [im (IRAC2) 
3.6 urn (IRAC1), 4.5 ftm (IRAC2) 
3.4 fim (Wl), 4.6 urn (W2) 
3.4 urn (Wl), 4.6 /_m (W2) 
3.4 /im (Wl), 4.6 fim (W2) 


23.6, 26.8 
23.6, 26.8 
100.lt 
177.lt 
84.7t 



tThe date quoted for the WISE images is the approximate median date of multiple exposures observed around that epoch. The Exposure time is simply the 

frame time (7.7 s) multiplied by the number of Lib frames selected for coaddition. 



TABLE C2 
Infrared photometry of SN 20()9js 



Dayst 


MJD 


Telescope 


Short MIR" 


Longer MlR b 








/uJy 


fUy 


-34.8 


55080.7 


Spitzer 


<9 


<18 


+2.2 


55117.6 


Spitzer 


226±13 


211±11 


+106.7 


55222.1 


WISE 


129±12 


136±15 


+295.4 


55410.8 


WISE 


<60 


<80 


+471.3 


55586.7 


WISE 


<90 


<100 



tDenotes days post-discovery (2009 Oct 11.44). 
"Refers to the Spitzer 3.6 fim or WISE 3.4 fim band. 
''Refers to the Spitzer 4.5 /im or WISE 4.6 ^m band. 
Quoted fluxes are directly-observed values, uncorrected for extinction. 
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